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SUMMARY 


A simple vortex system is used to model unsteady aerodynamic effects into 
the rigid-body longitudinal equations of motion of an aircraft. With the formu- 
lation used, only steady-state aerodynamic derivatives appear in the equations. 
It is found expedient to transform the equations into the frequency domain to 
make them useful for extracting aerodynamic parameters from flight data. The 
equations are used in the development of a parameter-extraction algorithm. If 
the algorithm is used with the unsteady aerodynamic modeling included, all 
extracted aerodynamic derivatives are the steady-state derivatives. If unsteady 
aerodynamic modeling is omitted, some extracted parameters will include the 
effects of unsteady aerodynamics and are interpreted as combinations of steady- 
state and acceleration parameters. Use of the two parameter-estimation modes, 
one including and the other omitting unsteady- aerodynamic modeling, provides a 
means of estimating some acceleration derivatives. Computer-generated data and 
flight data are used to demonstrate the use of the parameter-extraction 
algorithm. 


INTRODUCTION 

Extraction of aerodynamic parameters from flight data has become a routine 
task for numerous organizations. Sophisticated mathematical techniques have 
been developed for this purpose. Although parameter-extraction methods have 
been well developed and work very well with computer-generated data, several 
problems generally become apparent when applied to flight data. Some problems 
documented in reference 1 include: 

(1) Extraction of different numerical values for the same parameter, under 

similar flight conditions 

(2) Estimated variance in parameters which are much too optimistic when 

compared to values from ensemble averages 

(3) Dependence of values of extracted parameters on the shape of the 

control input 

One possibility that could have impact on these items is that unsteady 
aerodynamics associated with load buildup following an increase in angle of 
attack generally has not been included in parameter-estimation algorithms. 
Another possibility is that unsteady aerodynamics associated with downwash at 
the horizontal tail is generally approximated by assuming the existence of cer- 
tain acceleration derivatives, in particular the acceleration derivative C^. 
In the present study, the idea is that better consistency in extracted param- 
eters might be obtained if the two unsteady aerodynamic effects are modeled 
more precisely. 



Illllllllllffllllllllllllllilll I III 


Although various methods were available for calculating unsteady aero- 
dynamic loads on lifting surfaces (e.g., refs. 2 to 4), they were generally too 
complex for use in equations of motion or for parameter identification purposes. 
The study of reference 5 was a preliminary attempt to develop a method for 
including unsteady effects in equations of motion that could be used in param- 
eter identification. In reference 5 a simple vortex system was developed for 
estimating indicial lift and downwash associated with circulation for unswept 
wings in incompressible flow. This system gave results that were in good agree- 
ment with more accurate and complex vortex systems. The vortex system was 
generalized in reference 6, and the results were extended to obtain indicial 
lift for tapered, swept wings in incompressible flow. 

The purposes of the present study were as follows: (1) To develop further 

the methods of references 5 and 6 to derive equations for the downwash behind 
tapered, swept wings; (2) to develop equations of motion including unsteady 
aerodynamic effects; (3) to develop algorithms for parameter extraction account- 
ing for unsteady effects; and (4) to apply the algorithms to flight data. 
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SYMBOLS 

All aerodynamic coefficients are based on wing geometry, 
vertical acceleration, g units 
aspect ratio, b 2 /S 
elements of matrix 
span, m 
chord, m 

average chord, m 


lift coefficient. 


Lift 
qS 
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/ Cr , t D 




parameters defined in equations (28) , (29) , and (30) 


^m 


pitching-moment coefficient, ? A t c *1 i n gjom eiit , about aircraft center 

^ s w c w 

of gravity 


E ± ( ) 


f l' f 2' f 3 

F,G,H 


exponential integral of variable ( ) 


downwash factors 
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dv 


coefficients used in modeling downwash when indicial lift function 
is neglected 


acceleration due to gravity, 9.8 m/sec^ 
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ACj^ ( t ) 

Ae (t) 
e 
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moment of inertia about Y-axis, kg-m^ 

distance behind wing root quarter-chord point, m 

distance from aircraft center of gravity to aerodynamic center of 
horizontal tail, m 

aircraft mass, kg 

number of points 

pitch rate, deg/sec 

1 o o 

dynamic pressure, — pu S, N/m^ 

noise covariance matrix 
2 

area, m 
time, sec 
velocity, m/sec 
dummy variable 
downwash velocity 

distance downstream from aircraft center of gravity to wing aero- 
dynamic center, m 

constants in indicial lift function 
geometric angle of attack, deg or rad 
circulation strength, m 2 /sec 
deflection, deg or rad 
lift due to unit step in a 

response in downwash to unit step in a, deg or rad 
downwash angle, deg or rad 
parameter vector 
taper ratio, c t /c r 

sweep of quarter-chord line, positive for sweepback, deg or rad 
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air density, kg/m 3 
frequency , rad/sec 
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Subscripts : 
e elevator 

E effective 

f fuselage 

k variable index 

Z quantity at distance Z 

r root 

s unsteady aerodynamics only in downwash 

ss steady state (no unsteady aerodynamic effects) 

t horizontal tail or wing tip 

w wing 

Superscripts : 

variable in frequency domain 

* coefficient used in modeling downwash when indicial lift function 

is included 

• derivative with respect to time 
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Abbreviations : 


det 

determinant 

diag 

diagonal 

mag 

magnitude 

Re 

real part 


ANALYSIS 

The longitudinal aerodynamic parameters of an aircraft are determined pri- 
marily by the lift forces on the wing, horizontal tail, and elevator. In 
unsteady motion, the forces can be considered to be composed of three parts 
(ref. 7) . One part is due to circulatory flow associated with angle-of- attack 
variation where this variation is caused by the vertical (plunging) motion of 
the aerodynamic center. A second part is associated with circulatory flow and 
ascribed to the curvature of the flow streamlines relative to the surface when 
the surface is performing a pitching motion. This second component is often 
analyzed by considering an equivalent system which replaces the airfoil and 
circular streamlines with a circular-arc airfoil in rectilinear flow. The 
effective angle of attack associated with the curvature is the difference in 
surface slope between the quarter-chord and three-quarter-chord points. This 
effective angle of attack associated with rotation can be included with the 
angle of attack due to vertical motion for purposes of analysis (ref. 7) . A 
third part is due to instantaneous acceleration of noncirculatory potential 
flow. This lift is equal to the product of the virtual additional mass of the 
moving surface and the acceleration normal to the surface. This effect is 
generally small and is not considered in this paper. However, the effect could 
be important for very lightweight aircraft. 


Lift Due to Circulatory Flow 

The instantaneous lift on an airfoil associated with circulatory flow is 
obtained by use of an indicial lift function and the effective angle of attack. 
The indicial lift function AC L (t) is the lift response due to a unit step 
increment in the effective angle of attack. The analysis of reference 6 
resulted in a general expression for indicial lift associated with vertical 
motion, which could be fitted accurately by an exponential equation of the form 


AC ^ (t) = 

The constants y and z are functions of wing (or tail) geometry and are 
given in reference 6 for a wide range of wing geometries. The lift for arbi- 
trary variation in angle of attack can be obtained by using equation (1) in 
Duhamel 1 s integral as 


zut \ 
c r /2 j 

Ye / (1) 
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c L (t) 


(2) 


t 

Ac L (t - x)a E (x) dx 

0 

Application of equation (2) requires that the effective angle of attack be 
determined for a surface performing a plunging and/or pitching motion. The next 
two sections are concerned with determining the effective angle of attack for a 
wing and a horizontal tail surface. 


Effective Angle of Attack of Wing 

The flow around the wing and hence the wake produced by the wing are 
assumed to be uninfluenced by the presence of the horizontal tail. The effec- 
tive angle of attack of the wing can be considered to be associated with the 
flow normal to the wing at the quarter-chord point plus the curvature effect 
(mentioned previously) associated with pitching motion. The effective angle of 
attack can be shown to be 

(a E ) (3) 

' £i/ w u 2u 

where the last term is the curvature effect. The last two terms can be combined 
to yield 

K) w = a + u(* + ?) (4) 


Effective Angle of Attack of Horizontal Tail 

The effective angle of attack of the horizontal tail involves the same 
factors as that of the wing, but it is also influenced by the downwash from the 
wing. The effective angle of attack of the horizontal tail, therefore, is 

( a E ) t = a + ?( z t + t) " ( £ *) t (5) 

The approach used in determining downwash for arbitrary variation in angle of 
attack is to first calculate the downwash caused by a unit step increase in 
wing angle of attack, and then to use Duhamel's integral to obtain downwash for 
arbitrary variation in wing angle of attack. 

The simple vortex system illustrated in figure 1 and used in reference 6 
is used to calculate downwash. The system consists of a bound vortex along the 
wing quarter-chord line, a trailing vortex at each wing tip, and a shed vortex 
connecting the downstream ends of the trailing vortices. The vortices are all 
in a plane parallel to the free-stream direction, and the trailing vortices 
remain parallel to the free-stream direction. The vortex system stretches in 
the downstream direction at one-half of the free-stream velocity. A discussion 
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and justification for use of this vortex system to represent a complex physical 
situation is given in reference 6. The vortex strength is the same in all por- 
tions of the system. The time variation of the lift following a unit step 
increase in angle of attack and, consequently, the corresponding circulation 
strength are determined by equation (1) . 

In keeping with the simple vortex representation used, downwash is calcu- 
lated only in the plane of the vortex system, midway between the two trailing 
vortices. The downwash at a distance l downstream from the midpoint of the 
bound vortex can be determined by use of the Biot-Savart law (ref. 8) and is 
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In equations (7) , fj is associated with the quarter- chord- line bound vortex, 
f 2 (t) is associated with the two wing-tip vortices, and 13 (t) is associated 
with the shed vortex. In deriving equations (7) , use was made of the geometric 
relationship 

_b_ = A (1 + X) 
c r 2 

The Kutta-Joukowski equation used to relate lift to circulation is 


L(t) = pubHt) 

Solving for T (t) and substituting into equation ( 6 ) results in 


W l (t) = + f 2 (t) + f 3^ (8) 

The downwash angle is defined as 
W](t) 

e 2 (t) = — 

b u 


hence, equation ( 8 ) can be written as 


e Z (t) 


L(t) 

TTpu 2 bc r 



+ f 2 (t) 



Since the area of a wing is given by S = be = bc r 1 + the downwash equation 
can be written as 


e z<t> = + + f s (t 3 

The downwash following a unit step increase in angle of attack, therefore, is 

Ae z (t) = 1 4 + Tr X AC L (t) + f 2 (t) + f 3 (t)] (9) 


Substituting equation (1) into equation (9) results in 


A e z (t) 


1 + X 

4tt 



1 



[f x + f 2 (t) + f 3 (t)] 


( 10 ) 
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The vortex system being used is a simple model for a complex physical 
system; therefore, equation ( 10 ) is not expected to be very accurate quanti- 
tatively. However, the functional relationship between downwash and time 
should be useful. A means of improving the accuracy of equation (10) is to 
assure correctness of the equation at known conditions. In the case of down- 
wash, there are reasonably accurate methods of computing downwash under steady- 
state conditions. Under such conditions equation (10) becomes 

( Ae l)ss - + f 2 ( ”> + f 3<">] <H> 


Dividing equation (10) by equation (11) results in 


A e z (t) 




f i + f 2 (t) + f 3 (t) 
f-L + f 2 (~) + f 3 (») 


However , 



Therefore equation (12) becomes 



f x + f 2 (t) + f 3 (t) 
f l + f 2 (0O) + f 3 (°°) 


( 12 ) 


(13) 


The terms f 2 (°°) and f 3 (°°) can be evaluated by letting t -* °° in equa- 

tions (7) . The results are 


r 


l A ( 1 + A) 


A 


tan A 


f 2 <”> * I7TT 


TT ' 


l A (1 + A) 


tan A 


" mi + A) 

4 


+ 1 


> (14) 


f 3 (°°) = 0 


V 


Equation (13) is complicated because of the long expressions for f^(t), 
f 2 (t), f 3 (t), and f 2 (°°) as given by equations (7) and (14). After computing 

downwash for a few specific cases and some experimentation, it was found (as in 
ref. 5) that equation (13) could be approximated very closely by the equation 
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Equation (15) is also in a convenient form to use in frequency-response calcula- 
tions, as discussed subsequently, because Fourier transforms are available for 
the terms involved. Transforms are not known to exist for some terms of the 
downwash expression given as equation (10). Values of F 1 , G* , and H 1 are 
functions of wing geometry and length 2 and can be determined by curve fitting 
equation (15) to results calculated from equation (13) . Details of this pro- 
cedure are given in appendix A. 


Lift on Wing 

The lift associated with circulation for arbitrary variation in angle of 
attack of a wing (uninfluenced by the tail) is obtained by use of equations (1) , 
(2) , and (4) and is 



(16) 


Lift on Horizontal Tail 

The lift associated with circulatory flow of the horizontal tail is 
obtained by use of equations (1) , (2) , and (5) and is 



(17) 

The downwash for arbitrary angle of attack required in equation (17) is obtained 
by use of equation (15) in Duhamel's integral and is 
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Lift Due to Elevator Deflection 


In order to determine the lift associated with arbitrary variation of ele- 
vator deflection, an indicial response function must be obtained to use in 
Duhamel's integral. The lift would then be given by an equation of the form 


C Le <t) =y Ac L e (t " T)6 e (X) dT (19) 


where AC L ^(t) is the indicial function, that is, the time variation of lift 

associated with a unit step in elevator deflection. No data or convenient 
theory for obtaining ACL e were found to serve as a guide for developing a 
simple expression as would be required in parameter identification. However, 
on the basis of the wing indicial lift function, a control-surface indicial 
response might have the form 


AC L 


e 



( 20 ) 


The lift for arbitrary variation in control deflection would be given by using 
equation (20) in Duhamel's integral. 


EQUATIONS OF MOTION 
Unsteady Aerodynamics in Downwash 

I ndicia 1 lift function .- By assuming that initial pitch attitude and its 
variations are small, the rigid-body longitudinal perturbation equations of 
motion are 


PuS w 

a(t) = q ( t) C L (t) 


( 21 ) 


and 


q(t) = 


pu^S c 
K w w 

2I V 


Cm(t) 


( 22 ) 


The unsteady aerodynamic effects are modeled in the terms C L (t) and C m (t) , 
which are expressed as 


C L (t) = C L (t) + C L (t) + C L (t) 
w t e 


( 23 ) 
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and 


c " (t> ■ ■ + < v t> ] + ■ t CL » (t> 


(24) 


(All coefficients are based on wing geometry.) 

When equations (16) , (17) , (18) , and (19) are substituted into equa- 

tions (23) and (24) and the results are used in equations (21) and (22) , the 
consequence is a pair of integro-dif f erential equations involving convolution 
integrals. The solution of the equations is very time consuming even on a high- 
speed digital computer. Computation of sensitivity coefficients, which is an 
integral part of several parameter-extraction techniques, becomes very cumber- 
some. The equations, therefore, are transformed to the frequency domain to 
simplify them and make them more practical for use in parameter extraction. 

The equations of motion in the frequency domain can be shown to be 
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Equations (25) and (26) can be solved simultaneously to obtain the 
response characteristics as 


®°1<« ' 
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A il 

A 12 

-1 

puS 

w c 

2m L e 
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A 21 

A 22 


P^Vt ~ 

2I y L e 


where 


An _ 


PUS V 

2m 


K + ^0 


12 


= -1 + 


P s w 

2m 


x + 




+ 




2 


A 


21 


P u2s w 5 w 

2I y 


x 






f 



iOJ + 


P«S« 


21. 



+ 



xC T 


+ 



+ 



1 + c 


t^L + 


2 


( 29 ) 


(30) 


(31) 


f requency- 


(32) 


(33) 


(34) 


(35) 


(36) 
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In the parameter-extraction mode, flight data in the form of time histories 
of a, q, and < 5 e are converted to the frequency domain, and equations (25) 
and (26) are used in a parameter-extraction algorithm to extract parameters to 
best fit the frequency- domain data. However, some aspects of the equations are 
undesirable. First, the number of aerodynamic parameters in the equations is 


limited to 


( CL a) w ,ss' ( CLa )t,ss' (901)2,33' ( Cin a) 


and 


' L 6 f 


Some of these 


occur in combinations which multiply common states and, therefore, might be 
difficult to separate. Secondly, it would be preferable to extract aerodynamic 
parameters for the entire aircraft rather than those associated with aircraft 
components. These concerns led to the idea of examining the equations to deter- 
mine whether some alterations could be made to obtain parameters in a more 
desirable form. 


Unsteady aerodynamics are introduced into the equations of motion through 
the indicial lift and downwash equations. Results presented in reference 9 
indicate that effects of unsteady aerodynamics on aircraft motion could be 
modeled reasonably accurately if indicial response in lift is omitted and 
unsteady effects are included only in the downwash behind the wing. It was 
decided to proceed under this condition and to check at a later point to deter- 
mine whether the results justified making this simplification. 

Indicial lift omitted .- Including unsteady aerodynamics in the downwash 

and omitting the indicial lift function simplify equations (27) , (28) , (30) , 

and (31) . Note that F* , G*, and H*, which appear in equation (29) , also 

w w w 

involve the indicial lift (see eqs. (12) and (15)) and must be modified. Equa- 
tions (27) , (28) , (30) , and (31) become 


= ( c L a ) w , ss 

(SO. = ( CL *>t.ss 
( c h 2 ) s ■ ( c %) t , ss 

(%), ■ NeL 




(37) 

(38) 

(39) 

(40) 


where F , G , and H , which appear in the D-term, are obtained by curve 
fitting equation (15) to equation (12) with y = z = 0. (See appendix A.) 
Equations (37) to (40) are then used in equations (25) and (26) to replace 

C T , C T , C T , and C T , respectively. 

L w L t 1 h t 2 
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When unsteady effects are included only in the downwash f equations (25) 
and (26) with the use of equations (37) to (40) can be simplified and put into 
the forms 


and 



(41) 


(42) 


Some terms in equations (41) and (42) can be converted to more convenient forms. 
The first bracketed terms of equation (41) become, by using equation (38) , 



w, ss 


+ 




+ 



D 

ss 


(43) 


By assuming that the total lift of the airplane is contributed by the wing and 
tail 


(SJw.ss + ( C Ot, 


ss 


; - m 


l t ,ss 


= (%) 


(44) 


ss 


and therefore 


(%)„,„ + (%J S - K)ss + K) t , ss (H) v 


(45) 


ss 
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The second bracketed term of equation (41) can be multiplied and divided 
by c w /2 to ofc>tain 


iC + 1 K)w, 


+ - ^(Cl 

ss V C w c w /V 


'w)t,SS 


qc w = (c » 

2u l L q/ ss 2u 


(46) 


Some terms of equation (42) also can be combined. For example/ the first 
bracketed term can be combined with ( c m a ) f to obtain 


|;( Cl o)w.ss + Cfcti), + ( C %)f ‘ ( C "Vx)ss ■ £^( CL a)t,os(^) !f 


The second bracketed term of equation (42) can be shown to reduce to 


ss 


(47) 


c t\ Z t 


5 + ~jr;( CL «)w,ss + l z t + -rh-f* a) t . 


ss 


q _ 

u 


qc 


= r 

m q 2u 


(48) 


Using equations (45) and (46) in equation (41) , and equations (47) and (48) in 

equation (42), and the relationship — fC L g \ = /C mg \ results in 

c w ' e / S s ' e / ss 


and 


puS 


lOKX = q - 


w 


2m 


( CL a) ss + ( Clja )t, ss (sot ')i i ss ^ 


a 


c w 


Sl(%)ss 5 + ( CL «e) ss 6 « 


(49) 


icoq = 


P u2 S w c w 

21 v 


1 \ ~ 

( Ci %) 

ss 5 W ( CL a)t , ss \3a / 1 , ss 


a 


+ 25( Cm q) ss 5 + ( Cm 6 e ) ss <Se 


(50) 
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Equations (49) and (50) can be solved simultaneously to obtain the frequency- 
response characteristics as 



(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


As noted previously, the terms F w , G w , and H w are determined by curve 
fitting equation (15) to equation (12) with y = z = 0. A numerical procedure 
for this purpose and results for a wide range of wing geometries are discussed 
in appendix A. 

Equation (51) can be used for parameter extraction if the assumption that 
indicial lift effects can be ignored, as was assumed in the simplifications used 
to proceed from the complete equations (32) to (51), is valid. 
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Effects of Omitting Indicial Response From Equations of Motion 

The two steps taken to examine the effects of neglecting indicial response 
are as follows : 

(1) Determine whether the indicial lift has any effect on aircraft 
frequency-response characteristics . 

(2) Determine whether indicial response influences the numerical value of 
extracted parameters. 

The aircraft used in this part of the study is shown in figure 2 and has the 
mass and geometric characteristics of table I and the aerodynamic parameters of 
table II. This aircraft configuration was also used in flight tests discussed 
subsequently. The test aircraft had an all-movable tail for longitudinal con- 
trol; therefore for this configuration Cq^ = Cq^ , and the problem of deter- 
mining the indicial response for an elevator- type control did not have to be 
addressed. 

The computed frequency-response characteristics for the aircraft configura- 
tion with and without the indicial lift function, but including unsteadiness in 
downwash, are shown in figure 3. Since the frequency-response characteristics 
are very similar for the two cases, it appears that the indicial lift function 
would have a minor influence on aircraft motion, which is in agreement with 
the results of reference 9. In order to perform step (2), the two sets of 
frequency- response characteristics were used in a parameter-extraction exercise, 
using the algorithm developed in appendix B. The results are given in table III 
and show that the two sets of extracted parameters are in close agreement except 
for Cq . Since Cq has only a minor effect on airplane motion, the agreement 

in the remaining parameters seemed to justify further using the simplified equa- 
tion (51) , which neglects lift indicial responses, for parameter-extraction 
purposes. 


PARAMETER EXTRACTION 


Equation (51) , which includes unsteady aerodynamics only in the downwash, 
is used for parameter extraction. The procedure is to measure flight data as a 
function of time, convert the data to the frequency domain by using the discrete 
Fourier transform, and then use a parameter-extraction algorithm to estimate 
aerodynamic parameters to provide a fit to the data in some optimal manner. In 
the present study, the maximum likelihood algorithm outlined in appendix B was 
used for parameter extraction. In principle, it should be possible to extract 


the parameters 


and the product 


( C %) ss - K)ss' (Se)ss' ( C "«)ss' ( C ^) 


ss 


and 


ss 


( Cm <$e) 


SS 


( )t , ss^ a )~l f 


ss 
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Computer-Generated Data 


The frequency-response characteristics of figure 3, with unsteady effects 
only in the downwash, were used with two parameter-extraction algorithms: One 

included unsteady aerodynamics; the other neglected unsteady aerodynamics. In 
performing parameter estimation using the algorithm that included unsteady aero- 
dynamics, initial values of the parameters were assumed. These values were 
offset from the values which had been used to generate the data. The extraction 
program which included unsteady aerodynamics was allowed to iterate on the six 
parameters and the product previously listed. Although the extraction program 
retrieved the values used in generating the data, several parameters were very 
highly correlated. Pair-wise correlations above 0.95 were obtained between 

The probability of high correla- 

V ss 


( Crn a)ss' ( Cm q)ss' and ( CL a)t,ssGa) 


tions between these parameters had been anticipated because of the manner in 
which the terms appear in equations (49) and (50) . It had been hoped, however, 
that the frequency-dependent factor D would relieve the correlation problem. 
Since this did not occur, the extraction program was rerun with the product 


( CL a)t,ss(aa) Zt/ 


fixed at its correct value. This time the program again 


ss 


retrieved the correct aerodynamic parameters and all correlation coefficients 
were very low. It was therefore thought advisable to use this approach in sub- 
sequent parameter extractions to keep correlations reasonably low. 


The next step in this study was to use the same computer-generated data 
(unsteadiness only in the downwash) and the parameter-extraction algorithm which 
neglected unsteady effects (D = 0) . In this case it should be recognized that 
the unsteady aerodynamic effects would be absorbed in the extracted parameters. 
In particular the parameters approximated by 


and 




+ CL 




P S w c w C m d C L c 


4m 


(57) 


(58) 


would be extracted (ref. 1) . The second term on the right-hand side of each 
equation represents an approximation which accounts for frequency-dependent 
unsteady aerodynamics by use of constant factors. 

The parameters extracted with the two algorithms are shown in table IV. 
The two sets of aerodynamic parameters produce almost identical frequency- 
response characteristics as shown in figure 4. As indicated previously, the 
difference between /c m N and is interpreted as C^.. An approximate 

\ 1/ss 1 ^ 
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expression for c m^' i- s valid for low frequencies and is generally 

referred to as the lag-in-downwash effect (ref. 10) , is 



( 59 ) 


If the parameters of table II are substituted into equation (59) , the estimated 
value for is -7.33, which is very close to the difference between the 

extracted values of C m and C m shown in table IV (i.e., a difference of 

iiiq lu q 

-7.35). These results indicate that for this particular case unsteady aero- 
dynamic effects as approximated by the parameter C m . can be estimated accu- 

a 

rately by the lag- in- down wash effect. 

The reduced value extracted for when unsteady aerodynamics are not in 

the extraction algorithm is also of the expected order of magnitude. In this 
case the extracted value should be (appendix A of ref. 1) 


c ma ( Cm a)ss 


P s w c w c m^ C L a 

4m 


Using the geometric and mass characteristics of table I, the aerodynamic param- 
eters used in generating the data, and the value of C m . based on the geometric 

a 

parameters, C m ^ = -1.30, which is reasonably close to the value of -1.09 
actually extracted. 


Flight Data 

The flight data used in this part of the study were obtained for a general- 
aviation light airplane having the configuration shown in figure 2. The geo- 
metric mass characteristics and flight conditions are given in table I, and the 
measured flight data are shown in figure 5 as time histories of various measured 
states. Because of the measured vertical acceleration used in the parameter 
estimation procedure, the transformed equations of motion were extended by the 
equation 


a z = ^(icoa - q) 


As noted previously, it is much more convenient, when using a mathematical model 
involving unsteady aerodynamics, to perform parameter extraction in the fre- 
quency domain. The data of figure 5 were converted to the frequency domain by 
use of the discrete Fourier transform, and are presented as part of figures 6 
and 7. Parameter extraction was performed in two different modes: 
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(1) In the frequency domain, using the algorithm of appendix B with 
unsteady aerodynamic modeling included in the downwash. 


(2) In the frequency domain, using the algorithm of appendix B with 
unsteady aerodynamics not modeled- In this case some of the extracted param- 
eters, notably cl , reflect any unsteady effects. 


plus the product 


Mt,, 


ever. 


First, an attempt was made to extract all six parameters of equation (51) 

( ttM - The problem of high correlation was 

da h, ss 

encountered, as was the case with computer-generated data. In addition, how- 
some of the parameters were unrealistic in magnitude and/or sign. 

M)._ 

constant at an estimated value based on the wing and tail geometry. The param- 
eters extracted in these two modes are given in table V, and the frequency- 
response characteristics obtained by using each set of parameters are shown in 
figures 6 and 7. Both sets of parameters, when used in their respective model- 
ing equations, fit the flight-measured time histories equally well as judged by 
estimated residuals. 


eter extraction then was performed with the product 


( c %)t. 


Par am- 
held 


ss ' 


, ss 


Since it is preferable to extract individual derivatives, performing 
identification with modeling of unsteady aerodynamics as proposed herein appears 
to be preferable to extracting combination derivatives. Application of the 
extraction algorithm in the two modes suggested, that is, with and without 
modeling unsteady aerodynamics, provides a convenient means of separating the 
steady-state and acceleration derivatives. An alternate method of obtaining 
estimates of individual steady-state and acceleration derivatives by special 
flight maneuvers is suggested in reference 11. 


CONCLUDING REMARKS 

A simple vortex system has been used to model unsteady aerodynamic effects 
into the longitudinal equations of motion of an aircraft. The equations of 
motion in the time domain had two problems relative to application for parameter 
extraction: (1) Only a limited number of aerodynamic parameters appeared in the 

equations; and (2) the solution of integro-dif f erential equations led to lengthy 
computations. The second problem was circumvented by transforming to the fre- 
quency domain; however, the first problem (limited number of derivatives) 
remained. Subsequent calculations showed that unsteady aerodynamic effects 
were adequately accounted for when the unsteady effects were included only in 
the downwash. This condition permitted recasting the equations in the frequency 
domain so that the usual aerodynamic parameters were present in the equations. 

A parameter-extraction algorithm based on this formulation for unsteady 
aerodynamics permitted extraction of the aerodynamic derivatives associated with 
steady- state motion. If unsteady aerodynamic effects are omitted from the 
extraction algorithm, the extracted parameters will include the effects of 
unsteady aerodynamics; that is, some of the extracted parameters are inter- 
preted as combinations of steady-state and acceleration derivatives. Use of 
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the two parameter-extraction algorithms, one including and the other omitting 
unsteady aerodynamic effects, provides a means of estimating some acceleration 
derivatives. However, the acceleration derivatives obtained are constants that 
approximate the effects of frequency-dependent aerodynamic phenomena. 

The primary objectives of this study, the inclusion of unsteady aerody- 
namics in the rigid-body motion of an aircraft and the development of a 
parameter-extraction program which includes unsteady aerodynamic effects, have 
been accomplished. Additional studies suggested are: (1) To determine whether 

repeated tests under the same flight conditions will yield parameters with 
smaller ensemble variance when the extraction algorithm accounts for unsteady 
aerodynamics; and (2) to determine whether the extracted parameters will show 
less variation for different control inputs than is observed when unsteady 
effects are not taken into account. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 11, 1979 
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APPENDIX A 


CALCULATION OF F , G, AND H OF DOWNWASH EQUATION 

The downwash equation for a unit step increase in angle of attack for a 
wing represented by the vortex system of figure 1 is given by 




f ± + f 2 (t) + f 3 (t) 
f l + f 2 (oo) + f 3 (o °) 


(Al) 


The functions f 2 (t) , and f 3 (t) are defined by equations (7). As noted 

previously , it was found expedient to perform parameter identification in the 
frequency domain; however, this procedure led to another problem since equa- 
tion (Al) could not be converted because of the apparent lack of Laplace trans- 
forms for some of the terms involved in equation (Al) . It had been found, 
during the study of reference 4 , that for several representative wings, numeri- 
cal results from equation (Al) could be approximated very closely by an expres- 
sion of the form 




(A2) 


where F, G, and H are functions of wing geometry and distance behind the 
wing root quarter-chord point, Laplace transforms of the various terms of 
equation (A2) are readily available. 

In the present study, the constants F, G, and H were computed for a 
wide range of wing geometric characteristics by generating time histories of 
Ae 1 using equation (Al) , and then using a nonlinear least-squares procedure to 
determine values of F, G, and H which caused equation (A2) to provide the 
best fit to the time histories using the least-squares method. The procedure 
was to compute Ae ^ at 40 equally spaced intervals, starting at a low value 
of ut/c r and extending to a value for which Ae ^ was essentially constant 
(fig. 8 ) . These values were used as measured data y m . An initial guess was 
then made of a parameter vector 0 containing as its elements the constants F, 
G, and H. These initial guesses were used in equation (A2) to estimate Ae^/ 
which was treated as computed data y c . Sensitivity coefficients were computed 
analytically from equation (A2) and are 
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APPENDIX A 


A 


1 


_ 

3f 


l 


1 

1 - 





A 2 


3g 


Hut 

c r /2 




A 


3 


Hut 

9y c = r ut e c r/2 
3 h c r /2 


(A3) 


The sensitivity coefficients formed a row matrix A. An iteration procedure 
was then used to converge on the values of F, G, and H to best fit the 
measured data y m by using parameter updates computed from 

A0 = (2 A i A i) 

The updated parameter vector is 



0 = e Q + A9 (A5) 

where 0 q is the initial parameter vector. Iterations are continued until the 
parameter update is judged to be negligible or until the fit to the data meets 
some convergence criteria. This procedure was used to calculate the constants 
for a wide range of values of wing geometry (A, X, and A) and distance behind 
the wing root quarter-chord point. The results are given in figure 9. 


24 



APPENDIX B 


MAXIMUM LIKELIHOOD IDENTIFICATION ALGORITHM 

The parameter-estimation algorithm used in this study is the maximum like- 
lihood technique in the frequency domain (ref. 12). This method provides con- 
sistent estimates and has the asymptotic properties of unbiased and minimum 
variance estimates. The method which is based on maximization of the likelihood 
function yields the parameter-estimation algorithm and the covariance matrix for 
the parameters. In a separate step maximization of the likelihood function 
yields the covariance matrix for the measurement noise based on the current 
nominal solution. 

Let 0 be the vector of unknown parameters. The maximum likelihood esti- 
mation is based on maximizing the conditional probability density of z (n) 
given 0 or maximizing the log likelihood function L given as 


L (0 ) = N Re ^ £z(n) ~ x(n,0 o f] R 1 [z(n) - x(n,0 o )] 
n 

N i i 

- — log | R I - Constant 

= N Re ^ V(n) *R~ 1 V (n) - ~ log |r| - Constant (Bl) 

n 

where z (n) is the vector of measured variables and x(n,0Q) is the vector of 
computed output variables. 

To estimate R, the likelihood function is maximized with respect to the 
elements in R and thereby yields 


R = diag 



V(n) 


V*(n) 


(B2) 


For the estimates of the remaining parameters, the minimum of L(0) is found 
by using the modified Newton-Raphson method. The estimates are given by 


0 = 0q + A0 


(B3) 


where 0 q is the current estimate of parameter vector, A0 updates in the 
estimates of parameters, and A0 = -M“lg where M and g are, respectively, 
second and first gradients of the log-likelihood function. 
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APPENDIX B 


For fixed R = R 




v 


(B4) 


and 


3 2 L (0) 

M k l 96 k 36 z ~ RS 


2 


3v i 3v 
39^ R 30 


l 


(B5) 


That is. 


A0 = Re 


iZ o-i aiV 1 ! 

30 30 J 


V 3v* — l- 
Re Z 30- R V 


(B6) 


The effect of maximizing the likelihood function is the same as minimizing 
cost function (fit error) 


J = det 


iflisg 


V(n) V (n) 


Generally, in a convergent-estimation process, the fit improved with iteration 
as evidenced by a reduction in the cost function. Once the cost function 
settled so that change in two successive iterations, as defined by 
(J k “ Jk-lV J^, was less than 0.01, the parameters which maximize the likeli- 
hood function or minimize the fit error were considered identified. 

The maximum likelihood identification algorithm proceeds in the following 
manner : 

(1) Choose a nominal value for parameter vector 0q and calculate J(0q). 

(2) From the equations of motion obtain x(0Q,n) and the sensitivity 

3v (0Qf n) 

function — using finite-difference methods. 

(3) Compare the transformed flight data with the chosen nominal. 

(4) Calculate an update of 0 from the maximum likelihood estimate. 

(5) Calculate J(0) and compare with J ( 0 q ) . 

(6) Update the nominal parameter vector and continue until convergence. 
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TABLE I.- GEOMETRIC AND MASS CHARACTERISTICS OF AIRCRAFT USED TO 


STUDY VARIOUS EFFECTS OF UNSTEADY AERODYNAMICS AND USED TO 
GENERATE FREQUENCY- RESPONSE CURVES OF FIGURE 3 


Wing : 

Aspect ratio 7.35 

Taper ratio 1.00 

Sweep angle, deg 0 

Root chord, m 1.34 

Area, m 2 13.56 

Horizontal tail: 

Aspect ratio 4.21 

Taper ratio 1.00 

Sweep angle, deg 0 

Root chord, m 0.77 

Area, m 2 2.51 


Weight, N 
I Y / kg-m 2 
x, m . . 


1 1 , m 
1r ^Tl • * 

u, m/sec 
p, kg/m 3 


9230 

2135 

0.116 

1.34 

4.49 

4.38 

47.5 

1.076 
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TABLE III.- EFFECT OF INDICIAL LIFT FUNCTION ON EXTRACTED PARAMETERS 


USING COMPUTER- GENERATED DATA a 




Extracted parameter value 
(b) 


Parameter 

Data generated with 
unsteady downwash and 
indicial lift function 

Data generated with 
unsteady downwash but 
without indicial 
lift function 

° L a 

5.38 

(0.030) 

5.21 

(0) 

CL «e 

.67 

( -019) 

.74 

(0) 

% 

4.90 

( .600) 

11.02 

(0) 

Cm a 

-1.40 

( .006) 

-1.50 

(0) 

c m q 

-16.62 

( .250) 

-18.58 

(0) 

Cm 5 e 

-2.30 

( -014) 

-2.48 

(0) 


C . 33 


C . 33 



Parameter-extraction algorithm included unsteady aerodynamics in 
downwash but had no indicial lift function. 

^Values in parentheses are variances (Cramer-Rao lower bound) . 
c Value held constant during extraction. 
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TABLE IV.- RESULTS FROM PARAMETER EXTRACTION, USING COMPUTER-GENERATED 
DATA WITH UNSTEADINESS ONLY IN DOWNWASH 



a Values in parentheses are variances (Cramer-Rao lower bound) . 
^Value held constant during extraction. 
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TABLE V. - AERODYNAMIC PARAMETERS EXTRACTED FROM 


FLIGHT DATA OF FIGURE 5 


Extraction algorithm 
included unsteady 
effects in downwash 

(a) | 

Extraction algorithm 
neglected all 
unsteady effects 
(a) 

( CL a) ss 

5.12 (0.05) 

( Cl cc)ss = 

4. 

91 

(0.03) 

K)ss = 

4.69 (0.17) 

( CL q)ss = 

6. 

33 

(0.33) 

( CL O ss 

= 1.42 (0.07) 

( CLS *)ss 

= 1 

.14 

(0.05) 

( Cm a) ss 

-1.34 (0.01) 

C ™ol = • 

99 

(0. 

01) 

(c m _) 

-15.37 (0.01) 

= _ 20 

.77 

(0 

.58) 

V 3' SS 






( C "O ss 

= -3.25 (0.05) 

c m 6e - 2 

.95 

(0 

.01) 


(1^-) = b 0. 33 

\ da h, ss 


/ 9e 
\3a 

), 

= b 0 . 33 

SS 


a Values in parentheses are variances (Cramer-Rao lower bound) . 
*-Value held constant during extraction- 
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Figure 2.- Three-view drawing of test aircraft. 



++++++ Complete unsteadiness 
— Unsteadiness only in downwash 



Frequency, rad/ sec 
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(b) Pitch-rate response. 


Figure 3.- Concluded. 



++++++ Unsteadiness in downwash only 

Obtained using parameter extracted 

with steady-state model 
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(a) Angle-of-attack response. 


Figure 4.- Computer-generated data with unsteadiness only in 
downwash and data computed using parameters extracted by 
maximum likelihood method with unsteady effects neglected 
(steady-state model) . 
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sured from aircraft of figure 2. 


+ -*-•+•+■ + Flight data 
Computed 



Frequency, rad/ sec 

x io 2 



Figure 6.- Flight test data in frequency domain and data computed 
using extracted parameters when extraction algorithm included 
unsteady aerodynamics in downwash. 
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(c) Normal-acceleration response. 


Figure 6.- Concluded. 







++++++ Flight data 
Computed 




Frequency, rad/ sec 

(b) Pitch-rate response. 


Figure 7.- Continued. 
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Figure 9.- Continued. 
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